Full electrostatic control over polarized currents through spin-orbital Kondo effect 
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Numerical calculations indicate that by suitably controlling the individual gate voltages of a capacitively 
coupled parallel double quantum dot, with each quantum dot coupled to one of two independent non-magnetic 
channels, this system can be set into a spin-orbital Kondo state by applying a magnetic field. This Kondo regime, 
closely related to the SU(4) Kondo, flips spin from one to zero through cotunneling processes that generate 
almost totally spin-polarized currents with opposite spin orientation along the two channels. Moreover, by 
appropriately changing the gate voltages of both quantum dots, one can simultaneously flip the spin polarization 
of the currents in each channel. As a similar zero magnetic field Kondo effect has been recently observed by Y. 
Okazaki et al. [Phys. Rev. B 84, (R)161305 (2011)], we analyze a range of magnetic field values where this 
polarization effect seems robust, suggesting that the setup may be used as an efficient bipolar spin filter, which 
can generate electrostatically reversible spatially separated spin currents with opposite polarizations. 

PACS numbers: 72.15.Qm,72.25.-b,73.63.Kv,85.75.Hh 



Introduction. — Traditional spintronic devices rely on the 
use of ferromagnetic source and drain leads to produce and de- 
tect polarized spin-currents, like, for example, the Datta-Das 
spin field-effect-transistor fl]. More recently, the manipula- 
tion of single spins has become one of the paradigms for quan- 
tum information. To achieve easier integration with current 
technology, the use of semiconducting lateral single quantum 
dots (QD) has been suggested as a means to produce spin fil- 
tering and spin memory devices fl, which can be controlled 
through the use of electrostatic gates, without the need of fer- 
romagnetic contacts, nor highly inhomogeneous static mag- 
netic fields, or AC fields. Its experimental realization fl, us- 
ing a single QD and a large magnetic field to produce a bipolar 
electrically tunable spin filter, has spurred a multitude of pro- 
osals, e.g., two QDs embedded in an Aharonov-Bohm ring 
, a double QD (DQD) in parallel 15[], or in a T-shape geome- 
try |@], to cite a few. More related to the results presented here, 
Borda et al. J3l suggested the possibility of spin-filtering in a 
DQD device at quarter-filling, by exploiting spin and orbital 
degrees of freedom simultaneously through an SU(4) Kondo 
state. Right after that, Feinberg and Simon JH], by extending 
the ideas described in 0] to a similar DQD device, suggested 
the interesting possibility of a "Stern-Gerlach" spin filter ef- 
fect at half-filling. In this work we use two fully independent 
channels, and present detailed numerical results confirming 
the high efficiency of the spin filtering effect and suggest ex- 
perimental ways of observing it. 

The utilization of the Kondo effect jsjj] in a single-QD fioll 
has the potential to add an extra dimension to spintronics, as 
now the localized moment in a QD participates in a many- 
body state that may provide new functionalities to spintronic 
devices ll 111 . More complex Kondo-like regimes, like the so- 
called SU(4) Kondo state ITU [l3ll . may provide even addi- 
tional latitude to create, manipulate, and explore spintronic 
devices using QDs. In this work, we extend a recently ob- 
served variant of the SU(4) Kondo effect 11411 (dubbed the 
spin-orbital Kondo effect) to propose a device based on a ca- 
pacitively coupled parallel DQD which, when in the Kondo 



regime (through the application of a magnetic field — see 
below), functions as a bipolar spin filter that can produce 
currents with opposite polarities simultaneously (one in each 
channel of the DQD system). In addition, their polarities can 
be reversed by tuning the gate voltages of the QDs, i.e., the 
proposed bipolar spin filter is electrically tunable. As men- 
tioned above, a similar device had been suggested before Ji| . 
Here, we provide extensive numerical results to stimulate ex- 
perimental groups to try and observe this effect. 

Device and Hamiltonian. — The proposed setup is that 
of capacitively coupled parallel DQDs 01411 connected to 
completely independent metallic leads [see Fig. \Ha)] lfl5ll . 
Through an even-odd transformation, two leads decouple 
from the DQD and the system is reduced to that shown in 
Fig. 01b). Note that this transformation does not involve 
the QDs, therefore the interacting part of the Hamiltonian 
[(Eq. (2) below] remains unchanged. Then, the two-impurity 
Anderson Hamiltonian modeling our system is 
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The operator dl (d\,j) creates (destroys) an electron in QD 
A = 1, 2 with spin u = ±, while operator ct ■ (cAi+io-) does 
the same at site i in a non-interacting semi-infinite chain 
A = 1,2; n\ a = d\ a d\ a is the charge per spin at each QD, 
and both QDs have the same charging energy U. We include 
the effect of a magnetic field H acting just on the QDs fl. For 
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simplicity, we take the hybridization parameters t\ = t<i = 
t' . It is important to note that, contrary to Ref. JH, the only 
interaction between electrons in different channels A = 1 , 2 
is the inter-dot capacitive coupling V. Finally, as our setup 
consists of semiconducting lateral QDs, each of them can have 
different gate potentials V g \ and V g 2 [l5ll . and we will 
concentrate on the experimentally relevant regime U'/U < 
1.0 (the so-called SU{2) <g> SU{2) regime) iQ. All results 
shown were calculated using U as our unit of energy. The 
width of the one-body resonance for each QD is given by Y = 
nt' 2 po(Ep), where po(Ep) is the density of states of the leads 
at the Fermi energy Ep. Throughout the paper t = 1.0 and 
all the other parameter values are indicated in the figures or in 
the text. 



< b > u,v„, 





FIG. 1 : (Color online) (a) DQD connected to four metallic leads so 
that conductance through them can be measured independently. The 
QDs are subjected to an inter-(intra-)QD Coulomb repulsion U' (U). 
(b) After an even-odd transformation, two of the leads decouple and 
the system is reduced to just two leads coupled only through U' . 

This model has been studied extensively in previous works, 
and it is well known that for U'/U = 1.0 and zero-field it has 
an SU(4) Kondo fixed point 11211 . experimentally observed in 
Refs. M 1711 . Here, we want to address a completely differ- 
ent regime, although we will also show that our Density Ma- 
trix Renormalization Group (DMRG) 111 811 calculations faith- 
fully describe the SU(4) Kondo regime as well. 

Density Matrix Elements. — The results presented in this 
work were calculated using the DMRG fl9] and the Friedel 
Sum Rule (FSR) JHH^l. The validity of the FSR for the 
system studied here is discussed in the supplemental material 
023I1 . In order to characterize and identify different regimes, 
we use the reduced density matrix elements (DME), calcu- 
lated with the DMRG. The ground state wave-function can be 
written as 



l*o>=X>7,ffl7>l<*>. (5) 

7,(5 

where 7 stands for the 16 possible DQD configurations (0-0, 
cr-O, 0-cr, a-a' , 2-0, 0-2, er-2, 2-tr, and 2-2), while 6 represents 
the states associated with the Fermi sea. Summing over the 
band states 8 we obtain the weight projection of the different 
DQD configurations in the ground state. 
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As will be shown in Fig. [4] the diagonal DME can be used 
as a 'proxy order parameter' for the typical correlations that 



FIG. 2: (Color online) Conductance (in units of Go) for U' = 0.5, 
T = 0.02, H = 0.0 (a) and H = 0.05 (b). (m,n 2 ) specifies the 
occupancies of each QD [same values apply to (b)]. (b) H = 0.05 
suppresses the spin Kondo effect (note that color scales in each panel 
are different). We show that the bright (yellow) lines intercepted by 
the (white) dashed line, where G = Go, correspond to a peculiar 
Kondo effect. Points a and /? are discussed in detail in Figs. fj] and [4] 



characterize a many-body state like, for example, the Kondo 
state. This is very useful in the case of an unusual (or exotic) 
Kondo effect, where it may not be clear at first what are the 
relevant correlations that one should look for (from now on 
we generally refer to the diagonal matrix elements as DME 
weight, or simply DME). 

Numerical Results. — Figure 0a) shows the conductance 
G = G\ + G2 (G1.2 is the conductance for each channel) 
obtained through the FSR 12311 in the V 3 \ — V g 2 plane for 
U' = 0.5 at zero magnetic field. The different QD occu- 
pancies are indicated by the notation (m, 112). In Fig. |2jb) 
we present the conductance results for finite field H = 0.05, 
where the suppression of spin SU (2) Kondo in each channel 
can be clearly observed [color scales are not the same for pan- 
els (a) and (b)]. The (white) dashed line is the region of gate 
voltage variation in the V g i — V g 2 plane that interests us. It 
is parameterized by the expression V g 2 = —V g \ — (1 + 2U'). 
Conductance results along this line for 0.0 < H < 0.04 are 
shown in Fig. [5J where the (black) solid line shows results 
at zero magnetic field, with a well defined plateau around 
V g i = —1.0 [it corresponds to a cross section of the bright 
(yellow) region in Fig. 2(a)]. As the field increases, in steps of 
AH = 0.025 (dotted lines), the conductance at (and around) 
the particle-hole (p-h) symmetric point (V g i = —1.0) is sup- 
pressed very quickly, while narrow peaks start to form close 
to the charge degeneracy points [(red) dashed line], denoted 
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FIG. 3: (Color online) Effect of magnetic field over conductance 
along the (white) dashed line in Fig.^b). U' = 0.5, and V = 0.02, 
0.0 < H < 0.04 (field increases in steps of 0.0025). Half-filling 
conductance (around V g i ~ — 1.0) in the SU (2) ® SU (2) regime is 
suppressed faster than in the SU(4) regime (not shown). The dashed 
(red) line (H = 0.04) reaches unitary conductance G = Go at values 
of V 3 i corresponding to the points a and /3 in Fig. 2(b). 



a [(2,0)-(l,l)] and /3 [(0,2)-(l,l)] points in Fig. [2b), where 
G = Go . Note that these peaks are narrow along the (white) 
dashed line in Fig. [2jb), but along the charge degeneracy lines 
[the diagonal (yellow) bright lines in Fig. [2b)] they present a 
clear plateau structure. 

Now, in Fig. [4]we present one of the central results in this 
work. As an illustration of the use of DMEs to trace the possi- 
ble existence of a Kondo regime, panel (a) shows the DME 
for half-filling (ni + n-i = 2) configurations of the DQD 
for U' = 1.0 and zero magnetic field (SU(4) fixed point), 
as a function of V g i for V g 2 = —V g i — (1 + 2U') [equiva- 
lent to the dashed (white) line in Fig. |2b), but for U' = 1.0 
(see supplemental material) l23ll l. At the p-h symmetric point 
(V g i = V g2 = -U' - U/2 = -1.5) in Fig. Ha), i.e., at 
the half-filled SU(4) fixed point, one sees that the six possible 
two-electron configurations have all the same DME weight 
in the ground state, highlighting the fact that orbital and spin 
degrees of freedom are perfectly equivalent in the half-filled 
SU(4) Kondo state, i.e., spin and orbital degrees of freedom 
are maximally entangled I17II . This result is well known, but 
it serves to illustrate the use of the DME calculation to Took 
for' possible Kondo states. This is what is done in panel (b), 
where we present the DME results for U' = 0.5 and finite 
field H = 0.04 (same parameters as the ones for the (red) 
dashed line in Fig. [3]). In this case, we have two different val- 
ues of V g i for which we have two half-filling configurations 
with the same DME weight (same V g i values as the a and f3 
points in Fig . [3j ■ The crossing in the a (/?) point in panel (b) is 
between configurations t~t an d 2-0 (0-2). The important fact 
to note is that exactly at these crossings G = Gq (see Fig. [3j, 
indicating the possibility of a Kondo effect. 

Indeed, in Fig.[5]we show eight cotunneling processes (four 
in the upper panel and four in the lower) that shift the total S z 
spin of the DQD from S z = 1 to S z = (and vice-versa). 
The top processes correspond to the degenerate states f-t and 




FIG. 4: (Color online) (a) DME vs. VgiforV^ = -V g i-(1+2U'), 
T = 0.02, U' = 1.0, and H = 0.0. At the p-h symmetric point 
(Vgi/U = —1.5), we have a half-filling SU(4) Kondo regime, char- 
acterized by the entanglement of spin and 'orbital' degrees of free- 
dom, translated here into the equality of all the two-electron DME 
DQD configurations at V g \/U = —1.5. (b) Same as in (a), but now 
for U' = 0.5 (£77(2) <g> SU(2) regime) and H = 0.04, correspond- 
ing to the situation shown for the (red) dashed curve in Fig. [3] Note 
that the DME of the spin configuration f-j [(black) dotted curve] 
is the same as the 'orbital' configuration 2-0 [(red) solid curve] for 
V g i — —1.29 (a point), as well as for the configuration 0-2 [(blue) 
dashed curve] for V g i = —0.71 (j3 point). 



2-0 [a point in Fig. Efb)], while the bottom processes cor- 
respond to the degenerate states t~t an d 0-2 (j3 point). The 
virtual states contain either one or three electrons. The re- 
markable fact about these cotunneling processes is that they 
generate spin polarized currents in each channel, with oppo- 
site polarizations. In addition, once one sweeps V g i from a to 
(3, the polarization direction of the spin filtered current in each 
right-side lead is reversed: 4- (t) an d t (4-) m channels 1 and 
2, respectively, for the upper (lower) processes. Note that no 
other virtual states are connected (by t') to any of the degener- 
ate states [t-t and 2-0 (0-2)] in the a (/3) point J3. There are 
similarities between the Kondo effect described here and the 
one in Ref. Il4ll : from the lower inset on their Fig. 2 we see 
that the magnetic field raises (lowers) the energy of the config- 
uration l~l (t-f)> while maintaining the configurations t4--0, 
f-4-, and 4"t degenerate (to zero-order in t'). By adjusting 
the gate potentials V g i and V g 2, the configurations fi-0 and 
f-f can be made degenerate (a point). Then, the cotunneling 
processes in Fig. Ogive origin to the Kondo effect discussed 
here. 

In Fig.|6|a), we show conductance per spin type as a func- 
tion of Vgi (same parameters as Fig. [3] (red) dashed curve) 
for channels 1 and 2 (see legend). These results confirm that 



virtual virtual 




FIG. 5: (Color online) Top: Schematic representation of the four 
possible processes that flip the total spin from S z = 1 to S z = 
0, and back to 5 Z = 1, when the singlet state is in channel 1 (a 
point). The net effect is the transport of a spin down (up) in channel 
1 (2). The coherent superposition of these processes leads to Kondo- 
screening of the pseudospin associated to the two degenerate states 
with S z = 1 and S z — 0, and the consequent generation of a spin 
down (up) current in channel 1 (2). Bottom: Equivalent processes 
for the configuration where the singlet state is in channel 2 (J3 point). 
In this case, the polarity of the currents in channels 1 and 2 (when 
compared to top panel) is reversed. 
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FIG. 6: (Color online) (a) Channel conductance per spin vs V g i 
[same interval as in Figs.[3]and[4tb), for H = 0.04 and V — 0.2: 
at point a channel 1 (2) is polarized down (up), and at point j3 
channel 1 (2) is polarized up (down), (b) Polarization P\ for each 
channel calculated using data from (a), (c), (d) Polarization vs 
V g i (0.0025 < H < 0.04) for channels 1 and 2, respectively. 
Panel (c) also shows that for H/U = 0.025 there is a substantial 
range A V g i f» 0.045 for which the spin polarization varies between 
~ 70% and ~ 90%, indicating that the effect should be experimen- 
tally observable. 



the conductance at points a and ft are almost perfectly polar- 
ized, in accordance with the cotunneling processes described 
in Fig. [5] Figure 0b) shows the conductance polarization 
Px = (Gx t -G H )/(G xt + G xl ) for channels A = 1,2. 
Panels (c) and (d) show polarization results for channels 1 and 
2, respectively, for 0.0025 < H < 0.04 [other parameters 
as in panels (a) and (b)]. The results clearly indicate that the 
polarization effect is robust and does not require a high value 
of magnetic field. Indeed, as indicated in Fig.|3Jc), a polariza- 
tion of almost 90% can be achieved for H/U = 0.025. Taking 
U w l.OmeV for a GaAs QD|H results in H < LOT around 
the a point. As indicated by the double-head white arrow in 
Fig. He), a range of AV gl w 0.045, at H/U = 0.025, has a 
polarization varying from « 70% to « 90%. This indicates 
that there is enough range in the parameters space to allow for 
experimental observation of very high polarizations without 
the need of very high magnetic fields. 

Conclusions. — In summary, we have presented a peculiar 
Kondo effect involving a capacitively coupled parallel DQD, 
connected to two independent channels. To achieve this effect 
it is necessary to apply a moderate magnetic field and adjust 
the gate potential of each QD to take the DQD to a half-filling 
charge degeneracy point. The cotunneling processes in this 
Kondo effect are such that spin polarized currents are gener- 
ated in each channel, with opposite polarities. The analysis 
in Fig. indie ate that the effect should be experimentally ob- 
servable. 
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SUPPLEMENTAL MATERIAL 

Friedel Sum Rule and Conductance. — The first step in cal- 
culating the conductance through the DQD using the Friedel 
Sum Rule (FSR) |0, 2(1 21 ] is to use the following expression 
for the conductance [21] 



G\ a = —irTpxcr{E F ), 
h 



(7) 



where T is the coupling of each QD (A = 1, 2) to its cor- 
responding lead, and px er (Ep) is the local density of states 
(LDOS) of each QD at the Fermi energy. The FSR (valid if 
the ground state of the system is a Fermi liquid, which is true 
in our case as long as U' < U 1120 expresses a relationship 
between pxa{E F ) and n A(J (QD occupancy), 



P\a{E F ) 

where Ix is given by 

h = Im 



cos 2 (7r(n A(T -l/2) + J A ) 
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and G aa (cj) is the Green's function for each QD, 
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(8) 
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being E^ lb (w) the many -body self-energy, and E(cj) the one- 
body part of the self energy, which can be calculated exactly 
by taking U = 0. It can be shown 12211 that, for the parameters 
used in this work, the integral in the calculation of Ix is ap- 
proximately zero (within machine precision). Therefore, the 
conductance per channel, per spin, can be calculated as, 



G 



Act 



h 



sin 2 (7m A(T ), 



(11) 



and the total conductance, through both channels, is given by 



G — Gxa- 



(12) 



As an example, we will use eqs. (11) and (12) to calculate 
the conductance in the SU(4) regime (i.e., U' = U), by using 
DMRG to calculate the ground state charge per spin n A(J for 
each QD for a wide range of V g \/U values. We also apply a 
magnetic field that, contrary to previous works 111 31 . couples 
to the spins in the DQD, but not to the orbital degree of free- 
dom [see eq. (2)]. As there is complete symmetry between 
orbital and spin degrees of freedom in the SU(4) regime, we 
should obtain similar results to the ones obtained in Ref. fl 1 3f1 - 
Indeed, this is what is shown in Fig. |7Ja), where, at zero- 
field [(black) thick solid line], one sees the three characteristic 




FIG. 7: (Color online) (a) Conductance G (in units of 2e 2 /h) and (b) 
DQD total charge as a function of gate potential for U' — U (SU(4) 
regime), V g i = V g 2, and several different values of magnetic field 
H. In (a), the (black) thick line shows the conductance for H = 0, 
while thinner (black) dotted lines show the conductance for differ- 
ent H values in steps of -ff s tep — 0.0002. The thin (green) solid 
line shows the conductance for H = 0.0024, while the (red) thick 
dashed line corresponds to H = 0.05. For this large field value, the 
SU(4) Kondo effect at half-filling is totally suppressed and just the 
quarter-filled Kondo effect survives [total charge 1 and 3, see panel 
(b)]. These conductance calculations were done at T — through 
the Friedel Sum Rule [see eqs. (7) to (12)], by calculating the QD 
charges, as a function of gate potential, using the DMRG method. 
Note that a small T = 0.001 value was used to have the system well 
into the Kondo regime, thus resulting in very well defined charge 
plateaus and consequently very well defined conductance plateaus. 
Note that in panel (b) well defined plateaus are almost independent 
of the field H, although, as expected, in panel (c) large variations 
with field are seen, when the charge per spin is plotted as a function 
of gate potential. It is interesting to note that for H = 0.0024 the 
charge per spin for V g i/U — 0.0, —1.0, and —2.0 are similar (dips 
and peaks are indicated by arrows for V g /U = —1.0). This occurs 
because a small magnetic field is not enough to completely suppress 
the charge fluctuation (or mixed valence) effects. We noted that for 
a larger value of F the dips and peaks are even more pronounced, as 
expected (not shown), and for a much larger filed (H = 0.05, (red) 
dot-dash and open circle curves for spin up and down, respectively), 
they are completely suppressed. 



SU(4) conductance plateaus at quarter- and half-filling. For 
these calculations, we have used a small value of T = 0.001 
to place the DQD deep into the Kondo regime (at quarter- 
and half-filling). As shown previously II 1311 . the magnetic 
field strongly suppresses the Kondo effect at half-filling, but 
not at quarter-filling (see (black) dotted thin lines), in agree- 
ment with the fact that orbital and spin degrees of freedom 
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FIG. 8: (Color online) Spin [(red) solid line] and pseudo-spin [(blue) 
dashed line] correlations (see Ref. lfl2ll for definitions) compared 
with DME results as a function of U' /U, in panels (a) and (b), re- 
spectively. The parameters used were V — 0.02 and V g i — V g 2 — 
—0.5 — U' (i.e., p-h symmetric point). In panel (a), the SU(4) point 
(U'/U = 1.0) is where the curves for the correlations cross, while 
the corresponding point in panel (b) is where the DME for all the 
half-filling configurations are equal, reflecting the fact that spin and 
pseudo-spin degrees of freedom are perfectly entangled. Note also 
that the SU(4) point separates different regions: at left, one has an 
SU (2) ® SU (2) ground state, defined mainly by spin correlations, 
while at the right side one has a charge-ordered state, defined mainly 
by charge correlations fl2Tl . 



are equivalent in the SU(4) Kondo. The reason being that at 
quarter-filling one is left with an SU(2) orbital Kondo (in this 
case, where the magnetic field couples to spins, while one was 
left with an SU(2) spin Kondo in Ref. 11311 . where the field 
couples to the orbital degree of freedom). For H = 0.01 the 
SU(4) Kondo has been completely suppressed at half-filling 



[(red) dashed line]. Note that in panel (b) one sees very little 
field dependence in the total-charge variation with gate volt- 
age for the whole range of field explored (0.0 < H < 0.05). 
This is deceiving, though, as there is obviously considerable 
dependence with field in the gate voltage charge variation 
per spin, as shown in panel (c). In it, we see the interest- 
ing behavior at intermediate values of field (H = 0.0024, 
(green) dashed curve for spin up, and (green) solid dots curve 
for spin down), where, for gate voltage values for which the 
DQD passes through fluctuating valence regimes (i.e., close to 
charge degeneracy points), there are peaks (dips) in the spin 
down (up) occupancies (see arrows). Calculations (not shown) 
for larger values of T lead to more pronounced peaks and dips. 
Finally, note that for higher values of field (H = 0.05, for ex- 
ample), the peaks and dips have been suppressed [(red) dash- 
dot (open circles) line for spin up (down) in panel c]. 

Identifying regimes using DMRG Density Matrix Ele- 
ments. — As a way of demonstrating the use of the Density 
Matrix Elements (DME) to characterize different regimes in 
the DQD, we will compare DME results with those for spe- 
cific correlations that are usually used for that same purpose. 
It is known that at half-filling the DQD passes through three 
different fixed points as the ratio U'/U varies from zero to 
slightly above one: SU(2) <g> SU(2), for 0.0 < U'/U < 1.0, 
SU(4), for U'/U = 1.0, and a charge ordered state for 
U'/U > 1.0 1 1 211 - Figure [S^a) shows the evolution of spin 
[(red) dashed line] and pseudo-spin [(blue) solid line] corre- 
lations as U'/U varies from 0.0 to 1.6 (see Ref. 1 1 211 for def- 
initions). First, concentrating on the SU(4) point, one sees 
that both correlations are identical at that point, which is in- 
dicated by a dashed vertical line that separates the two re- 
gions associated to the other two fixed points. For small val- 
ues of U'/U the spin correlations dominate, the pseudo-spin 
correlations being very small. The opposite occurs for the 
U'/U > 1.0 region, indicating the charge ordered state 1112 1. 
In Fig. [8jb), we show the corresponding DME results for all 
the half-filling configurations. The results neatly agree with 
the overall picture obtained through the correlations in panel 
(a), demonstrating once more the use of DMEs to uncover dif- 
ferent phases in strongly correlated electron models. 



